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1  Introduction 

This  manuscript  discusses  the  recent  results  in  the  research  project  on  inviscid  regularization 
techniques  with  emphasis  on  applications  to  compressible  gas  flows.  It  is  well-established 
that  general  systems  of  quasilinear,  first-order  hyperbolic  PDEs  exhibit  finite-time  blowup 
of  smooth  solutions,  thus  restricting  wcll-posedness  of  classical  solutions  in  short-time.  For 
instance,  with  conservation  laws,  one  must  generally  seek  global-in-time  solutions  in  a  larger 
class  of  discontinuous  solutions.  Further  physically  motivated  admissibility  conditions  are 
placed  in  order  to  obtain  the  uniqueness  and  stability  of  such  solutions. 

The  goal  of  such  techniques  is  to  modify  the  equations  by  smoothing  out  the  nonlinear 
terms  through  a  mollification  or  convolution  process.  Consequently,  higher  regularity  to 
the  solutions  is  achieved  while,  at  the  same  time,  capturing  the  physical  behavior  of  the 
original  model.  Mathematically,  however,  such  qualitative  properties  must  be  verified  in 
order  for  this  technique  to  be  accepted.  More  specifically,  when  applied  to  a  system  of 
partial  differential  equations  usually  supplemented  with  appropriate  initial  data,  the  global 
wcll-posedness  of  smooth  solutions  should  be  shown.  Moreover,  the  limit  of  the  solutions  for 
the  filtered  system  should  converge  in  some  sense  to  a  physical  or  entropy  solution  of  the 
original  system.  In  the  subsequent  sections,  we  provide  the  class  of  functions  that  will  be 
used  as  the  spatial  filters  in  our  method  along  with  the  analytical  results  we  have  obtained 
when  applying  filtering  to  both  the  transport  model  and  for  the  more  general  quasilinear, 
symmetric  hyperbolic  systems  of  PDEs. 

2  Averaging  kernels  and  Filters 

Let  us  describe  the  appropriate  Liters  we  shall  implement  in  the  averaging.  For  a  given 
real-valued  function  /,  we  introduce  an  averaging  kernel  G  and  define  the  filtering  of  /  by 
the  convolution  operator 

J(x)  =  G*f(x)=  [  G(x-y)f(y)dy 
Jn 

on  the  given  domain  R  C  MV 

Given  some  integer  k  >  1,  we  only  consider  averaging  kernels  G  G  IUfc_1’1(lRAr)  that  satisfy 
the  following: 
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Properties 

Mathematical  Expression 

Normalized 

JmNGdx=  1 

Non-negative 

G(x)  >  0  for  all  x  6 

Symmetric 

|xi|  =  |x2|  =>  |G(xi)|  =  |G(x2)| 

Non-increasing 

|xi|  <  x2|  =►  |G(xi)j  >  G(x2)| 

Table  1:  The  properties  of  the  averaging  kernels 


In  the  physical  sense,  the  averaging  should  only  provide  non-negative  weight  to  particles, 
have  no  preferential  direction,  and  should  give  more  weight  to  particles  that  are  physically 
closer.  In  addition,  we  prescribe  a  parameter,  a  >  0,  to  such  a  filter  such  that 

Ga  =  —G  (- 

a  \a 


This  parameter  a  acts  as  a  scaling  of  the  kernel  and  controls  the  level  of  filtering.  One 
example  of  a  commonly  studied  filter  is  the  Helmholtz  filter  /  =  /  —  a2fxx  corresponding  to 
the  averaging  kernel 


Ga(x)  =  —  exp  -J-f 


2a 


\x\ 


a 


3  Transport  equations 

Consider  the  multi-dimensional  transport  (or  transportation)  equations  which  adjoins  the 
continuity  equation  with  the  inviscid  Burgers  equation 

r  pt  +  v  •  (pu)  =  o 

<  ut  +  u  ■  Vu  =  0  (1) 

l  (p(x,0),u(a;,0))  =  (p0(x),u0(x)). 

The  transport  equations  provide  a  simple  system  of  two  conservation  laws  for  which  classical 
weak  solutions  may  not  exist  for  general  initial  data,  thus  a  distribution  solution  consisting 
of  Dirac  delta  functions  (delta-shock  solutions)  must  be  introduced  [1,2].  On  the  other  hand, 
the  transport  equations  model  the  dynamics  of  particles  that  adhere  to  one  another  upon 
collision  and  has  been  studied  as  a  simple  cosmological  model  for  describing  the  nonlinear 
formation  of  large-scale  structures  in  the  Universe  [3].  For  example,  u  may  represent  the 
flow  field  carrying  dust  particles  with  density  p,  and  the  delta-shock  wave  represents  a 
concentration  of  dust  on  a  shock  which  attract  the  dust  [4].  The  filtering  method  employed 
is  the  observable  divergence  method  and  the  initial  value  problem  (1)  is  modified  as  follows 

{pt  +  pV  ■  u  +  u  ■  Vp  =  0 

ut  +  u\/u  —  (2) 

(p(x,  0),  u(x,  0))  =  (po(x),  u0(a;)). 

The  following  existence  result  is  shown. 
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Theorem  3.1.  Let  p0  :  fl  ^  M.  and  uq  :  12  — s ►  M>N  be  bounded  C 1  functions,  then  there  exist 
unique  classical  solutions  u(x,t )  and  p(x,t )  to  the  IVP  (2). 


4  Quasilinear  symmetric  hyperbolic  systems 

In  this  work  we  proposed  a  regularization  method  for  quasilinear  symmetric  hyperbolic  sys¬ 
tems  through  a  filtering  of  the  coefficients  and  the  source  terms-similar  to  the  filtering 
described  above.  The  primary  motivation  here  is  to  generalize  and  unify  the  previous  re¬ 
sults  on  the  Burgers’  and  homentropic  Euler  equations  while  encouraging  the  application  of 
similar  filtering  techniques  to  other  physical  models  other  than  those  found  in  the  study  of 
compressible  gas  dynamics. 

Let  us  consider  a  symmetric  hyperbolic  system  of  N  equations  in  n-space  variables 


n 


(3) 


where  the  kb’s  are  symmetric  N  x  N  matrices  while  h,u  are  iV-vector- valued  functions.  We 
always  prescribe  an  initial  condition  to  this  system 


u(x,  0)  =  u0(x). 


(4) 


Throughout,  we  also  assume  that  the  system  is  strictly  hyperbolic,  i.e.  each  N  x  N  matrix 
Ai(x,t,u )  has  N  distinct  eigenvalues 


A«  <  A<‘>  <  ...  <  A«. 


We  introduce  spatial  averaging  to  the  coefficient  matrices  in  (3)  to  prevent  the  finite¬ 
time  blowup  of  solutions,  thus  providing  the  global  well-posedness  of  smooth  solutions  to 
our  modified  IVP  (provided  we  have  sufficiently  smooth  initial  data).  More  precisely,  given 
a  filter,  our  spatial  averaging,  regularization  technique  modifies  (3)  into  the  following: 


n 


u(x,  0)  =  u0(x)  for  i6K” 


(6) 


where  the  bar  represents  the  convolution  product  taken  with  respect  to  the  filter  G  in  the 
^-variables.  More  precisely, 
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The  averaging  for  the  vector-valued  function  h  is  defined  similarly.  Moreover  we  impose  the 
condition: 

max  (||M(x,f,M(a;,f))!|i:0O,  ||/i(a;,£,u)||ij00)  <  const  (7) 

for  any  solution  u  G  C([0,T],  Hk  n  C1((0,  T),  Hk~x)  for  (5)-(6).  where  u  =  G  *  u  is  the 
convolution  product  taken  with  respect  to  the  x-variables.  First,  we  briefly  compile  several 
standard  results  on  linear,  symmetric  hyperbolic  systems  required  in  proving  the  main  exis¬ 
tence  theorems  for  the  quasilinear  case.  The  reader  is  referred  to  [5-7]  for  further  details  on 
their  proofs. 

Linear  symmetric  hyperbolic  systems 

The  linear  system  on  Euclidean  spaces  has  the  form 


vt 


+  ^2  A(z,  t)vXi  =  h(x,  t )  in  UT 


(8) 


i= 1 


with  initial  value  u(x,0)  =  u(0)  =  v0(x).  The  following  conditions,  we  refer  to  by  (Cl),  will 
be  imposed  on  the  linear  system. 

(a)  k  >  1  +  n / 2  and  vq  G  Hk. 

(i b )  Ai  are  symmetric  . 

(c)  t  ^  Ai{t)  =  M-,t)  is  of  class  C([0,T],Hk(Rn,RN2)). 

(d)  t  ^  h{t)  =  hi-,  t)  is  of  class  C([0,  T],  Hk(Rn,  RN)). 

Proposition  4.1.  Suppose  that 

veCi^TlH^nC1^),^) 

satisfies  the  initial  value  problem  to  (8),  then  v  satisfies  the  energy  estimate 


max  ||u(f)||n  <  e 

o  <t<T 


CT 


1^(0)111  +  2  /  11^)112* 


(9) 


where  the  constant  C  depends  on  the  supremum  of  A  and  its  first-order  spatial  derivatives 
on  Ut- 

Another  more  general  estimate  is  the  following. 

Proposition  4.2.  Suppose  that 

v  G  C([0,  T},Hk)  n  ^((0,  T),  H k~k) 

satisfies  the  initial  value  problem  to  (8),  then  u  satisfies  the  energy  estimate 


max  ( 

o  <t<T 


k(*)IU  +  IhMIU-i)  <  cte&T  ( ||ii(o)||t  +  / 


(10) 


where  the  constants  C\;  and  [3k  depend  on  the  supremum  of  A  and  its  spatial  derivatives  up 
to  order  k  on  Ut- 
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Proposition  4.3.  The  initial  value  problem  to  the  linear  system  (8)  has  a  unique  solution 
of  class  C([0,T},Hk)nC\(0,T),Hk-1). 


Quasilinear  symmetric  hyperbolic  systems 

For  the  quasilinear  case  the  following  conditions  will  be  made  and  will  be  referred  to  by 
(C2).  Define  Br  C  Hk  to  be  the  closed  ball  with  radius  R. 

(a)  k  >  1  +  nj 2  and  uq  G  Hk. 

(■ b )  For  u  G  Hk ,  A(x,t,u)  and  h(x,t,u )  are  //^-functions  that  satisfy  (Cl). 

(c)  The  maps  u  G  BR  i— >  Ai(x,  t,  u )  and  u  G  BR  i— ►  h(x,  t,  u ) 

are  bounded  (maps  bounded  sets  to  bounded  sets)  and  are  C 1  maps  with 
bounded  derivatives. 

Theorem  4.4.  Suppose  we  have  the  a  priori  bound:  Given  any  solution  u  G  C1(MTl  x  [0,  T]) 
of  the  IVP  (3)-(4),  there  exists  a  constant  CR  >  0  depending  only  on  T  >  0  such  that 

sup  {\DaA(x,t,u(x,t))\,  \Dah(x,t,u)\}  <  Ct-  (11) 

UT,\a\<k 

Then  the  IVP  has  a  unique  classical  solution  in  C'1(Mn  x  [0, T]). 

The  proof  is  given  in  three  main  steps.  In  step  1,  we  set  up  an  approximate  iteration  of 
linear  systems  along  with  a  corresponding  transformation  related  to  the  global  solutions  to 
these  linear  systems.  In  step  2  we  show  that  this  transformation  is  a  strict  contraction  on 
an  appropriate  function  space  for  sufficiently  small  time.  Further,  the  unique  fixed  point  of 
this  contractive  map  agrees  with  the  unique  short-time  classical  solution  to  the  quasilinear 
IVP.  Step  3  applies  energy  estimates  along  with  (11)  to  extend  this  unique  classical  solution 
further  in  time. 

Proof.  (Theorem  4.4)  Step  1:  Choose  an  arbitrary  T  >  0.  We  shall  prove  existence  of 
solutions  up  to  this  arbitrary  time.  First  construct  the  linear  problem: 

vt  +  Ai(x ,  t,  u(x,  t))vx.  =  h(x,  t,  u)  in  UT  (12) 

v(x,  0)  =  u0(x).  (13) 

where  the  subscript  i  is  short-hand  for  summation  from  1  to  n.  The  global  existence  and 
uniqueness  of  solutions  v  G  C([0,  T\,  Hk )  D  C((0,  T),  Hk~1)  to  this  IVP  holds.  The  first  step 
to  showing  existence  of  a  solution  to  the  quasilinear  system  is  to  consider  the  transformation 
T  defined  by  v  =  Tu  where  u  is  given  and  v  is  the  solution  to  (12)-(13).  Our  goal  is  to  prove 
this  transformation  is  a  strict  contraction  on  a  suitable  function  space.  We  consider 

u  G  Xk’T  =  C([0,  r],  Hk{Rn,  RN)). 
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Using  the  energy  estimates,  one  has 

max  (\\v(t)\\k  +  IM^IU-r)  <  K1eIC2T  ^||v(0)||fe  +  ||h(s)||fcch^  ,  (14) 

where  the  constants  K,  depend  only  on  the  constant  Ct  from  (11).  Dehne 

B =  {ue  Xk’T  :  ||u||xfclT  <  R}. 

It  is  clear  from  (14)  that  T  maps  to  itself  for  sufficiently  small  r  and  a  suitable  R.  We 
now  show  that  T  is  a  contraction  on  B^T  in  the  X0,T-norm. 

Let  Vj  =  Tuj  for  j  =  1,  2  and  set  w  =  v\  —  v2.  Then  w  satishes  the  linear  system 

Wt  +  Ai(u\)wx  =  H(x,t)  and  w(0)  =  0 

where  H(x,  t )  =  h(x,  t,  u i)  —  h(x,  t,  u2)  +  {A^x,  t,  u2)  —  A^(x,  t,  Ui))(y2)x.  From  the  Lipschitz 
continuity  of  h  with  respect  to  u  and  Sobolev  embedding,  ||//’(i)||0  <  C\\u\  —  xt2 1| o  where 
the  constant  C  depends  on  R  and  the  Lipschitz  constants  of  Ai  and  h.  Using  the  energy 
estimate  (9),  we  obtain 

max  \\Tui  —  Tu2\\i  <  CeK2Tr  max  llwi  —  rt2 II n ■ 

o<t<r "  u  o<t<r  u 


Hence  T  :  Bk^T  i— >  Bk^T  is  a  strict  contraction  with  respect  to  the  X°’T-norm  for  sufficiently 
small  r. 

Consider  the  iteration  scheme:  let  u^+1)  =  'TvP'1  with  =  uq.  As  a  consequence  of 
the  contraction  mapping  principle,  tb-7')  converges  to  a  unique  u  G  X0,r  i.e. 

lim  max  II —  u|L  =  0.  (15) 

j—>00  0<t<T  ' 


Step  2  :  We  show  in  this  step  that  this  limiting  function  u  belongs  in  C1(MTl  x  [0,  r],Kw). 
Energy  estimates  and  interpolation  inequalities  imply  that,  for  any  s  with  0  <  s  <  k, 


max  \\uU)  — 

0  <t<T 


<  Ck  max  || 

0  <t<T 


u  llo 


max 

0<t<r 


—  U 


(l)\\s/k 

life 


<  C  max  || 

0  <t<T 


-  7/U)  II  1~s/fc 
u  llo 


(16) 

(17) 


It  follows  from  this  and  (15)  that 


lim  max  ||«^  —  u\\s  =  0 

j— >00  0<t<T 

for  any  0  <  s  <  k.  Choosing  s  so  that  s  >  1  +  n/2,  Sobolev  embedding  implies 

w(i)  C([0,  t],  C1(Mn)). 
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In  addition,  it  follows  immediately  from 

(u^  —  u^-%  =  h{xR,u^^)  — 

+  A(x ,  t,  u ^  ~2^)ujf  -1)  —  A{x ,  t ,  -u^  _1'))m^7 

that  u ^  converges  to  u  in  C1([0,  r],  C'(Mn)),  thus  u  G  C1(Mn  x  [0, 7~] ,  M^) . 

Step  3:  In  this  step,  we  extend  the  local  classical  solution  to  the  whole  interval  [0, T\. 
The  condition  (11)  and  the  previous  a  priori  energy  estimate  imply  that  we  have  a  uniform 
bound  max0<t<T  ||u(t)||fc  <  At  for  some  constant  KT  >  0.  Therefore,  we  can  take  the 
radius  of  the  ball  R  to  be  sufficiently  greater  than  Kt  then  repeat  the  above  local  existence 
argument  on  (r,  2r),  (2r,  3r),  (3r,  4r), . . .  until  we  have  covered  [0,  T],  □ 


Application  to  spatially  averaged  systems 

Let  us  apply  theorem  4.4  to  spatially  averaged  systems.  To  do  so,  we  verify  condition  (11) 
holds.  Given  |a|  <  k ,  Young’s  inequality  and  (7)  imply 

\\Daa^\\00  <  ||G||fc-iIi||aJ])||i)0O  <  \\G\\k-hlconst. 

A  similar  argument  holds  for  h ,  thus  verifying  (11).  Thus  we  have  shown  the  following. 

Corollary  4.5.  Consider  the  IVP  to  the  regularized  system  (5)-(6)  along  with  the  same 
assumptions  made  on  Ai}  h,  and  uq  as  was  stated  in  theorem  4-4  and  condition  (7).  Then 
there  exists  a  unique  global-in-time  classical  solution  for  this  IVP. 

5  Zero  a  convergence  to  weak  solutions  conservation 
laws 

Our  inviscid  regularization  technique  has  an  important  application  to  conservation  laws.  We 
want  to  study  the  convergence  of  the  solutions  ua  as  a  limits  to  zero  in  order  to  justify  such 
a  regularization  method  in  the  sense  that  it  captures  the  ‘behavior’  of  the  original  equations. 
For  instance,  several  important  questions  arise.  Does  the  sequence  of  averaged  solutions 
converge  in  some  appropriate  function  space  as  the  filtering  parameter  tends  to  zero?  If  so, 
does  this  limit  function  satisfy  the  original,  non-averaged  system  in  the  weak  sense?  Consider 
the  following  Id  system  of  conservation  laws  with  forcing 

ut  +  f(u)x  =  h(x,  t,  u )  in  f/T  =  lx  (0,  T)  (18) 

u(x,  0)  =  u0(x)  for  (19) 

where  u ,  /,  h  are  M>N  vector-valued  functions  on  Ut ■  Let  A  =  Df  and  place  the  same 
conditions  on  the  resulting  quasilinear  system 

ut  +  [Df}(u)ux  =  h(x ,  t,  u )  (20) 
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as  was  given  in  theorem  4.4.  In  addition  to  these  previous  assumptions,  we  place  further 
conditions  as  follows. 


Further  Assumptions:  Let  n  —  1,  k  —  2,  supj/  \h(x,t,0)\  <  C0  for  some  constant  C0,  the 
initial  data  has  bounded  total  variation,  the  Jacobian  of  the  flux  /  is  a  diagonal  iV-by-IV 
matrix,  A(u )  =  diag(Xi(u),  X-2 (u), ...,  Xn(u)),  and  let  the  Liter  satisfy 

ll<(=“-““)IU  ;„  =  °(£)  (2D 


for  all  u  E  L\oc 

We  apply  our  inviscid  regularization  to  this  system  of  conservation  laws  to  obtain 

ut  +  A(ua)ux  =  ha(x,t,ua).  (22) 

However,  we  consider  a  slightly  different  type  of  averaging 

ut  +  A(ua)ux  =  ha(x,t,ua).  (23) 

Observe  that  if  one  has  flux  function  /  =  f(u )  that  is  a  quadratic  polynomial  in  u  (e.g. 
Burgers’  and  homentropic  Euler  equations),  then  A  is  linear  in  u  and  the  two  averaged 
systems  are  equivalent.  Fortunately,  for  general  smooth  flux,  the  global  existence  theory  for 
the  Id  system  (23)  still  holds. 

Our  aim  for  this  section  is  to  verify  that  the  global  classical  solutions  for  the  Cauchy 
problem  to  (23)  does  in  fact  converge  to  a  weak  solution  for  the  Cauchy  problem  to  (18)- 
(19).  By  a  weak  solution  we  mean  a  solution  to  the  Cauchy  problem  in  the  following  sense. 

Definition  5.1.  A  function  u:lx[0,r]H  M>N  is  a  weak  solution  of  the  Cauchy  problem 
(18)-(19)  if  u  is  continuous  as  a  function  from  [0,  T]  into  L]oc,  the  initial  condition  (19) 
holds  and  the  restriction  of  u  to  the  open  strip  Ut  is  a  distributional  solution  i.e. 

poo  poo 

I  u<pt  +  f(u)<f>x  +  h(xU,u)(j)  dxdt  +  /  u0(x)(f)(x,  0)  dx  =  0  (24) 

— oo  J  —  oo 

for  every  C°°  function  (f  with  compact  support  contained  in  the  set  M  x  (— oo,  T). 

The  notion  of  proving  the  convergence  result  is  summarized  in  two  key  steps.  In  step 
1,  the  required  uniform,  BV,  and  L 1  estimates  are  established  on  the  sequence  of  averaged 
solutions  {F"}a>o  that  guarantee  compactness  in  C'([0,  oo),  L(oc(M.)).  In  step  2  the  limit 
function  in  the  a  — >  0  limit  is  shown  to  satisfy  the  definition  of  a  weak  solution  for  the 
Cauchy  problem.  Without  loss  of  generality,  we  prove  the  case  for  the  scalar  case  N  =  1. 
The  proof  follows  similarly  for  N  >  1  as  a  consequence  of  the  diagonal  assumption  we 
imposed  on  the  matrix  A. 

The  following  highlights  the  results  needed  in  obtaining  our  desired  convergence  to  weak 
solutions  result. 


Proposition  5.2.  The  sequence  of  solutions  m“:Kx[0,T]h1  to  the  Cauchy  problem  to 
(23)  satisfy  the  following  bounds. 

| ua(x,  t)  |  <  Mi  for  all  x,  t , 


T.V.(ua(-,t ))  <  M2  for  all  t 
where  Mi  and  M2  are  independent  of  a. 


With  these  a  priori  bounds,  we  can  then  show  the  following. 

Theorem  5.3.  The  corresponding  sequence  of  averaged  solutions  ua  :  M  x  [0,  T]  i— >  M. 
from  Proposition  5.2  has  a  subsequence  {w7}7  that  converges  strongly  to  some  function  u 
in  C'([0,  T\,  Lloc(M))  as  7  limits  to  zero.  Moreover,  this  limit  function  u  is  a  weak  solution 
to  (18)-(19). 


We  should  remark  that  depending  on  the  filter  in  question,  further  conditions  on  the 
flux  /  may  be  necessary  in  order  for  (21)  to  hold.  To  illustrate  this,  consider  the  Helmholtz 
filter.  As  the  subsequent  calculations  will  show,  we  must  take  /  =  f{u )  to  be  a  quadratic 
polynomial  in  u. 

Using  the  definition  of  the  Helmholtz  filter,  we  may  multiply  (23)  by  a  test  function  0 
then  integrate  by  parts  to  obtain 


ua<ft  +  f(ua)4>x  +  ha(x,t,ua)(J)dx  dt  +  /  uo(x)(p(x,  0)  dx 

Jr 

f*T  /*  2  pT  /* 

=  /  /  f'(ua) (u°  -  ua)(t)xdxdt  -  %-  /  /  f"(ua)(uf.)2(f)xdxdt 

Jo  Jr  2  J o  7m 

-  W  r  /  r^)(KT<t>dxdt 

z  Jo  Jr 
=  Ei  +  E2  +  U3. 


Clearly,  E\  limits  to  zero  as  a  tends  to  zero  since  \\ua  —  ua\\Li  =  a2||u“a,||ii  =  0{1  /a)  for 
the  Helmholtz  filter  and  E3  =  0  by  our  restriction  on  /.  It  remains  to  be  shown  that  the 
term  E2  limits  to  zero  in  the  a  zero  limit.  This  follows  from  the  following  estimate 

V  [T  !  f"(.ua)K)2<f>*dxdt 
z  Jo  Jr 

a2  /  u"  dx  dt 

Jo  Jr 

<  Q||/,,(Ma)IU»||^lk<x.||Ma|UccM2T^  a. 

Again,  compactness  of  the  sequence  of  solutions  will  allow  us  to  formally  take  the  limit 
as  a  tends  to  zero  and  the  limiting  function  is  a  global  weak  solution  to  the  original  Cauchy 
problem. 


—  ^  ll///('uQ)IUool|0a;||L°o  ||w“||z,o 
2  V  ct 
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6  Numerical  Simulation 


Our  regularization  technique  has  been  quite  successful  in  regularizing  Burgers  equations  and 
homenetropic  Euler  equations  [8-14],  In  the  following  we  present  simulation  results  for  our 
ID  regularized  Equation. 

ID  Regularized  Euler  Equations.  The  equations  for  the  Euler  equations  without  the 
homentropic  assumption  are 


Pt  T  P^x  T  ^ Px  0 

(25a) 

(pu)t  +  puux  +  u(pu)x  +  Px  =  0 

(25b) 

(pe)t  +  peux  +  u(pe)x  +  Pux  +  uPx  =  0 

(25c) 

P  =  (7  -  1)  (pe  ~  \pu2^\  • 

(25d) 

With  analytical  techniques  these  equations  have  been  proven  to  converge  to  weak  solutions 
of  the  original  equation  provided  modest  assumptions  on  the  solutions.  We  have  also  found 
analytical  traveling  wave  solutions  for  both  sets  of  equations  which  has  led  to  analytical 
results  showing  that  shock  thickness  can  be  controlled  by  varying  the  amount  of  filtering 
used  in  the  low-pass  filter. 

Numerically  the  equations  are  shown  promising  behavior.  The  solution  appear  to  be  well 
regularized  and  capture  much  of  the  behavior  of  the  original  equations.  Figure  1  shows  a 
double  shock  tube  simulation  for  the  modified  homentropic  Euler  equations  plotted  against 
the  solution  to  the  original  homentropic  Euler  equations.  Both  the  expansion  wave  and 
shock  behavior  are  being  captured.  Figure  2  shows  a  shock  tube  simulation  for  the  Euler 
equations.  The  modified  equations  are  showing  a  regularized  solution  that  is  capturing  the 
expansion  wave,  contact  surface,  and  shock. 
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Figure  1:  This  figure  shows  a  numerical  simulation  of  the  modified  homentropic  Euler  equa¬ 
tions  (dashed  line)  plotted  against  the  solution  to  the  homentropic  Euler  equations  (solid 
line).  Here  the  value  of  a  =  0.05.  In  both  figures  it  is  clear  that  the  modified  homentropic 
Euler  equations  are  capturing  both  the  expansion  wave  and  shock  behavior,  (a)  The  velocity, 
(b)  The  density. 


Figure  2:  This  figure  shows  a  numerical  simulation  of  the  modified  Euler  equations  (dashed 
line)  plotted  against  the  solution  to  the  Euler  equations  (solid  line).  Here  the  value  of 
a  =  0.05.  In  the  figures  it  is  clear  that  the  modified  Euler  equations  are  capturing  both  the 
expansion  wave,  contact  surface,  and  shock  behavior,  (a)  The  density,  (b)  The  velocity,  (c) 
The  energy. 
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